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ABSTRACT 

Tidal and tidal-resonant effects in coalescing compact binary systems are in¬ 
vestigated by direct numerical integration of the equations of motion. For the 
stars polytropic models are used. The tidal effects are found to be dominated 
by the (non-resonant) /-modes. The effect of the g-mode-tidal resonances is 
obtained. The tidal interaction is shown to be of interest especially for low- 
mass binaries. There exists a characteristic final plunge orbit beyond which 
the system cannot remain stable even if radiation reaction is not taken into 
account; in agreement with results obtained by Lai et al. ( |1993 ). The impor¬ 
tance of the investigated effects for the observation of gravitational waves on 
Earth is discussed. 
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1 INTRODUCTION 

Close neutron star binaries are among the primary sources of gravitational waves that will be 
detected by the future gravitational-wave observatories on Earth ( [Thorne 1987|) . Using the 
matched hltering technique proposed by Thorne (tl987|) , close binaries should be detectable 
out to cosmological distances in the Gpc regime. The observed gravitational waves will 
be produced during the hnal inspiral and coalescence stages of the motion of the both 
compact objects (neutron stars and/or black holes). The detection process will last 2 - 
3 secs as the gravitational waves are expected to be measured when the binaries will be 
in the gravitational-wave-frequency window of 100 - 1000 Hz. Future improved detector 
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2 Kostas D. Kokkotas and Gerhard Schafer 

technology is expected to widen this window down to 10 Hz and up to 10 kHz, which means 
longer observation time and higher resolution of the signal, leading in this way to a more 
accurate determination of the system parameters like the masses, spins, and radii of the 
bodies. Nevertheless, already the 100 - 1000 Hz window can give us important information 
about the above parameters of the system if we know the orbital motion and the interaction 
mechanism of the two extended bodies with reasonable accuracy (during the final inspiral 
the bodies cannot be considered as point-like). The better our theoretical knowledge will be 
the more information will be extractable from the detected signal while any effect that will 
be omitted will reduce the signal to noise ratio of the detector. For example, the omission of 
post-Newtonian orbital terms in the construction of the template waveform will not change 
significantly its amplitude but it will produce destructive effects on the phase of the wave 
leading to a decrease in the correlation of the incoming wave and the applied template 
( ^chutz 1993|) . Recent numerical simulations for testing the post-Newtonian approximation 
( [Kokkotas et al. 1994a|) have shown that the correlation falls down by ~ 40% for nearly 
equal-mass systems while the phase effects become more dramatic for systems with unequal 
masses. Moreover, this phase mismatch will produce analogous effects in the accuracy of the 
calculation of the system parameters (e.g. masses) as well as of the observational distance 
and, as a consequence, of the rate of the detected events ([Kokkotas et al. 1994a]) . For these 
reasons a lot of attention has been paid recently to the study of possible fluid characteristics 
that may corrupt the signal ([Kochanek 1992]; [Jaranowski fc Krolak 199^ ; [Bildsten fc CutlCT 
1992]; ISpyrou fc Kokkotas 1994[) . 


What is expected from theoreticians in the remaining 5 to 6 years before operating 
gravitational-wave observatories might be available is the qualitative and quantitative study 
of every aspect of the compact binary-star coalescence in order to construct the best possible 
theoretical models that our theoretical approximations allow ( Cutler et al. 19"^ ; Schutz 
1993[) . Tidal effects have been proposed ([Fabian et al. 1975[) as the reason for the formation 


of close binary systems. Tidal capture has been extensively studied by various authors in 
the past two decades leading to a quite well understanding of the influence of the tides not 
only on the formation but also on the evolution of close binary systems ( Press fc leukols'ky 
[1977] ; [Lee fc Ostriker 198^ ; [Alexander 1987[) . 

In this line of research we are investigating tidal and tidal-resonant effects in the final 
stages of the inspiral of two neutron stars. Beyond this we investigate the information that 
can be gained about these effects from the detected gravitational waveforms. We also calcu- 
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Tidal and Tidal-Resonant Effects in Coalescing Binaries 3 
late the evolution of the binaries in the Ist-post-Newtonian (IPN) approximation to general 
relativity, including 3.5PN gravitational radiation damping, and compare it with the tidal 
evolution. 

In the next section we present equations that govern the evolution of close binary systems 
under the influence of Newtonian tidal and 2.5PN gravitational (quadrupole-)radiation re¬ 
action forces. In the third section we describe numerical procedures and discuss the obtained 
results. 


2 EQUATIONS OF MOTION 


Let us consider a close binary system consisting of two bulk-non-rotating neutron stars with 
masses mi and m 2 . Both of the stars are more or less tidally distorted by their respective 
companions. The Hamiltonian governing the motion of the two extended bodies, in our 
case, consists of four parts: the contribution from the stationary orbital motion, a part 
describing the oscillations of the stars, a term arising from the gravitational interaction of 
the tidally deformed star 2 with the undeformed star 1, and vice-versa, and, hnally, a part 
which describes the gravitational radiation reaction or radiation damping. 

The orbital Hamiltonian, at Newtonian order, in polar coordinates reads 


1 


Horb = — \Pf+^ 


Pl\ G^M 


2/i 


i?2 


R 


( 1 ) 


where R denotes the separation of the two stars and <I> the orbital rotation angle. Pr and 
are the corresponding canonical momenta. 
mim2 


T = 


and M = mi -|- m 2 


( 2 ) 


mi -|- m2 

are the reduced mass and the total mass of the system, respectively. 

The Hamiltonian describing the oscillations of the stars corresponds to an inhnite set of 
harmonic oscillators. It has the form ([Alexander 1987|): 




( 3 ) 


where are generalized coordinates and Pn their conjugate momenta; a;„ are real frequency 
eigenvalues. The index a on Hose denotes the stars; it takes the values 1 and 2. Like in eq. 
(3) we often drop indices for simplicity. 

The tidal gravitational potential at a point (r, 9, 0) inside, let us say, the star 2 (the 
centre of mass of this star is put into the origin of the coordinate system) due to the star 1, 
which is located at {R, 7r/2, d>) and treated as point-like, reads in spherical coordinates: 
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AirGmi “ ^ 

^T = -— 2 ^ 2 ^ 


R f^2 r^-l 2^+1 


YU'^l2,^)Yfje,f>){rlR)\ 


( 4 ) 


where are normalized spherical harmonics. Using this expression, the gravitational po¬ 
tential energy \Jt of the deformed star 2 (with the corresponding mass change dp) in the 
tidal held of star 1 is given by 

= J (frdpd^r. (5) 

By the use of the continuity equation (po is the unperturbed mass density) 

dp = -divpo^ (6) 

and integration by parts, the expression (5) becomes 

= j Poi- Vct>Td\ (7) 

where ^ is a vector describing the Lagrangian displacement of the huid. The vector can be 
written as a linear superposition of normalized eigenvectors 

= I]grx(t)^„(r), (8) 

n 

which satisfy a self-adjoint eigenvalue equation of the type 

= Po^nin- ( 9 ) 

The normal modes are characterized by their spherical harmonic indices ^ and m, and 
the radial index n. They can be decomposed into radial and poloidal spherical harmonics 
([Chandrasekhar 1964|): 


in = </>) = + Llir)rV\ 17 „^( 0 , (j)). 

By the aid of the unitary transformation 


ph) _ _ p it 

^n£m ' ^nE—r 


1 


j(o) 

•^nim 


^ni—r 


(0 < m < i), 
{1 < m < i) 


( 10 ) 

( 11 ) 

( 12 ) 


iv^ 

we can create an orthonormal real basis in which the generalized coordinates {a = e, o) 
are real. Herewith Ut can be written as 


= - 


oo e G) , -.i 

n 0 O rr—P r> nrv^ —H '2 ' ' 


Gm 17712 


n 1=2 (T=e,o m=0 
rO = pn«lrnd>'l rricrOld)! = 


(13) 


where trig(^^(<h) = cos(77i<h), trig^^($) = sin(m<h). <h decomposes into $ = n -|- p, where v is 
the true orbital anomaly and g = constant. The prime on the fourth sum means exclusion 
of the case m = 0, a = o, and 72 is the radius of the star 2. Qnt denotes the overlap integral 
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pr 


t+i 


dr, 


s 


(14) 


where r = r/r 2 and p = pQr\/m 2 is a dimensionless density and Wirn is given by 




= {-)- 


dvr 


2 £ + 1 


(£ — m)!(£ + m)! 


' i -\- m 


m 


(15) 


where (—)^ is zero if k is not an integer. The fnnctions qlf2n of motions of the harmonic 


oscillator, dehned by the Hamiltonian (3), in terms of action-angle variables, and 
respectively, can be represented as 




(16) 


where 0^^^ = ojnpt + (fnem — constant. Herewith, hnally, f/^takes the form 




X 


Gl7liTTl2 yjj 1 f^2 

r) /? O /T=P n vy, —n '2 ' 


Qni^ 


- 1/2 

ni 


n i= 2 ( 7 =e,om =0 
(JnLy^Hcosiel^!^ + m$) + cos(-0iE + 


+ sin(-^SL + 


3 ( 0 ) 


(o) 


(17) 


Circnlar orbits play an important role in coalescing binaries. There, <h is given by the 
single Keplerian orbital angnlar freqnency ujk- The orbital motion indnces pertnrbations on 
the star modes which rotate with the phase velocities mux, see eq. (13). For an adiabatic 
inspiral sitnation ujk is a slowly varying fnnction with time. Whenever ojni/oJK approaches 
integer valnes, resonances can occnr; see eq. (17) and, for remember the condition 

(£ + m)/2 = integer. 

The Hamiltonian of the entire system - orbital motion, star oscillations, tidal conpling 
of the two stars, radiation reaction - takes the form: 

H[t) = + Hill + hSI + Hj]f' + (18) 


where — U^’^\ ete. 

The (explicitly) time-dependent reaction Hamiltonian Hreac 
I99§: 


is given by, cf. (|S chafer 


Hreac (f') 


2G d^Qijjt) 

dC 


P^Pj 


GMp 


~W 


(19) 


where Qij = p {jTR^ — denotes the mass-qnadrnpole tensor of the two-body system. 

A spin-orbit-interaction part Hso conld have been also inclnded into the Hamiltonian, 
e.g. see ( Pamonr fc Schafer 1988|) . Bnt it has been fonnd that in onr case of bodies withont 
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6 Kostas D. Kokkotas and Gerhard Schafer 


bulk rotation it does not contribute significantly to the evolution of the system, and in what 
follows we shall not discuss it any further. 

The system evolves in time as dictated by the Hamiltonian equations of motion: 


dP 

dH{t) 

dH 

dH{t) 

dt 

aR ’ 

dt 

dP ' 

dpn 

dH{t) 

dqn _ 

1 dH{t) 

dt 

dqn 

dt 

dpn 


( 20 ) 

( 21 ) 


The respective reaction parts, Fp 
take the form, cf. ( [Schafer 1985| ): 

Fpn 

Fp, = - 


Fr = 


= 


SG^Pr f 

GM^u 

3 [ 

5 

8 / 

'2GM^u^ 

5 (FR^v \ 

R 

8 

(2P^ + 6 

15 FR^u 

SG^PrP^ 

3 FuR^ 




and Fr^, in the above equations, in polar coordinates. 



( 22 ) 

(23) 

(24) 
(26) 


The gravitational waves are emitted primarily by the orbital motion. The amplitudes of 
the leading (quadrupole) waves are proportional to 1/i? and their phases vary with time 
through 2<h. 


3 RESULTS AND DISCUSSION 


The Hamiltonian equations of motion (20) - (21) have been integrated for various binary sys¬ 
tems consisting of neutron stars from the less relativistic to the extreme relativistic regime. 
The equation of state was taken to be of polytropic type. 


p = Kp^, where T = 1 -|- 


N 


(26) 


with polytropic index N = 0.5,1,2 while the central density was ~ 10^“^ — lO^^gr/cm^. It 
is known that the best approximation for constructing realistic polytropic neutron stars 
is to take polytropes with indices 0.5 up to 2 ( [Finn 1987]) . The adiabatic index Ti we 
have chosen in such a way that gf-modes exist. This is possible as far as the adiabatic 
index Ti is larger than T; in the limit Ti = T there are no g-modes. There are two other 
important aspects regarding the selection of the appropriate adiabatic index: As Ti —> T 
the g-mode frequencies are slowing down significantly and the g-modes become degenerate 
at zero frequency. This property of the g-modes in principle increases the influence of the 
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Tidal and Tidal-Resonant Effects in Coalescing Binaries 7 

interaction part of the Hamiltonian (since Ut ~ but there are also the values of the 

overlap integrals which decrease much faster as it can be seen from Tables 2a and 2b. The 
result is that the interaction Hamiltonian becomes less important as Ti —T. For a detailed 
discussion of the g-modes we refer to Finn ( |1987|) , McDermott et ah (|1985|) , and Reisenegger 


and Goldreich (|1992|) and references therein. For the /-modes it does practically not matter 
which values of Fi we use, the frequency and the overlap integral are not very sensitive 
to the choice of Fi. In what follows we shall discuss polytropic stars which show typical 
external characteristics, i.e. masses and radii of accepted relativistic neutron stars models. 
More general equations of state will surely alter the results, but the main outcome of the 
paper, in particular those effects which are dominated by the /-modes, should have wider 
applicability. We have decided to study the binary evolution for nine different star models, 
A, B, C with polytropic indices 0.5, 1, and 2 (the value of the constant K has been chosen 
appropriately in each case). These models span the whole range of the neutron star family, 
from the very condensed state with high surface redshift to the less compact state with small 
surface redshift, see Table 1. In all cases the equations of structure are Newtonian-like, i.e. p 
in eq. (26) is the Newtonian mass density. Thus, with the models we have constructed, the 
less relativistic neutron stars are the better approximated ones. For these models normal¬ 
mode eigenfrequencies and eigenfunctions have been found and overlap integrals have been 
calculated. For i = 2, the part of the spectrum around the /-mode together with the hrst few 
gravity and pressure modes for the N = 0.5,1, 2 polytropes are shown in Tables 2a and 2b 
where additionally the corresponding values of the overlap integrals are listed. Although for 
our numerical calculations we have mainly used the values of Table 2a we present also Table 
2b in order to show the dependence of the g-mode frequencies and the operlaping integrals on 
the values of the adiabatic index. The frequencies agree with those given by Cox ( p.98CI|) . For 
an easier use of Tables 2a and 2b, in the last column of Table 1, the fundamental frequencies 
are given in units of Hz. When we checked the procedure for the calculation of the overlap 
integrals we have calculated also the normal modes and the overlap integrals for N = 3 
polytropes (Fi = 5/3). The agreement with the results of Lee & Ostriker (|1986|) was very 
good, usually the difference was less than 0.1% even for the lowest ^f-modes for which special 
care must be taken in order to avoid numerical errors. 

The procedure followed for the calculation of the normal-mode eigenfrequencies is similar 
to Dziembowski (|1971|) with the exception that the hnal frequencies and eigenfunctions 


© 0000 RAS, MNRAS 000 , 000-000 














8 Kostas D. Kokkotas and Gerhard Schafer 

calculated by the aid of this procedure have been used as trial expressions for a variational 
calculation of the eigenfrequencies. Practically, this final step together with an extra check 
for the orthogonality of the eigenfunctions was just a verification for the accuracy of our 
method. 

The Hamiltonian equations of motion (20) - (21), for each value of n, i, m {m positive), 
and cr, form a system of 6 first-order ODEs. As can be seen from the form of Ut (see eq. (17)), 
because of the existence of the term (r 2 /i?)^, respectively (ri/i?)^, only the harmonics with 
7 = 2 and 3 are important (m = 0,..., 7). Furthermore, from the Tables 2a and 2b follows 
that the /-mode together with the first g- and p-modes maximize the overlap integral, but 
since Ut is inverse proportional to the square root of eOni the p-modes do not contribute 
significantly and thus their incorporation in the following calculations is not needed. 

The initial conditions put two non-oscillating stars into circular orbits with orbital fre¬ 
quencies of b^Hz (correspondingly, the frequencies of the emitted gravitational waves are 
\QQHz). Notice that at the late stages of the evolution of a binary system orbits are expected 
to be circular because of the tendency of the gravitational radiation damping to circularize 
the orbits. Then the system is left to evolve according to the equations (20) - (21). The 
evolution of the system depends on the radiation damping and on the tidal interaction of 
the two stars. 

Since the important aspect of this work is the comparison of the rate of evolution of 
binary systems with and without tidal interaction we have to discuss also a refinement of 
the binary motion where the components are treated as point-like. For this we consider the 
IPN orbital dynamics. In case of circular motion the (secular) radiation-reaction change of 
the IPN orbital angular frequency reads, e.g see ([Cutler fc Flanagan 1994|) 


duj _ 
dt 





(27) 


In this relation the gravitational radiation damping has been taken into account up to the 


3.5PN level (for non-circular orbits, see ([Junker fc Schafer 1992|) ) using balance equations. 


A local expression for the 3.5PN reaction force is still not available. 

As seen from Figure 1, the evolution of the system as calculated by the damped IPN 


approximation (equation 27) is slower than the evolution which assumes damped Newto¬ 


nian approximation. If we additionally turn on the tidal part of the Hamiltonian, Had, then 
the evolution becomes faster since the orbital motion will excite oscillations of the initially 
non-oscillating stars. The amount of the energy transferred from the orbital dynamics to the 


© 0000 RAS, MNRAS 000 , 000-000 










Tidal and Tidal-Resonant Effects in Coalescing Binaries 9 
stellar dynamics depends on the stellar eqnation of state and the excitation efficiency of the 
tidal conpling. As mentioned earlier, the overlap integral Qn 2 is maximnm for the /-mode, 
and it has been shown numerically that all other modes together cannot absorb more orbital 
energy than the /-mode alone (see also ( [Reisenegger 199^ )). If radiation reaction is turned 
off, the tidal interaction makes the orbits oscillating around the initially chosen circular or¬ 
bit. Without radiation damping the system is a near-integrable system since the part of the 
Hamiltonian which describes the orbital motion and the oscillations of the star is integrable 
while the other part which describes the tidal interaction is the “perturbation Hamiltonian”. 
Therefore, techniques for the treatment of evolution and resonances of near-integrable sys¬ 
tems can be applied ([Alexander 1987| ; [Goldstein 1980| ; [Lichtenberg fc Lieberman 1983[) . For 
an extensive analysis of the resonances in Newtonian binary systems the reader is referred 
to a paper by Alexander ([1987[) . 

For the cases we are studying, i.e. coalescing binary systems consisting of two neutron 
stars, strong /-mode resonances can occur only when the two stars are near their hnal plunge 
stage. Since the orbit there shrinks very fast, the passage through the tidal resonances is 
instantaneous and thus it does not result in a signihcant energy exchange between the quasi- 
periodic orbital part and the periodic internal part of the system. But if we allow for the lower 
frequency ^f-modes then resonances can occur as it is shown in Figures 2 and 3. Obviously, 
for ujn 2 = i-e. £ = 2 and m = 2, the tidal-interaction Hamiltonian leads to strongest 
resonance, see eq. (17), and the ( 7 -mode frequencies are low enough to fulhl this condition, 
see Table 2. However, the Figures 2 and 3 show that the /-mode (non-resonant) effects are 
much more important than ( 7 -modes resonant ones. The effect of the resonances depend on 
the adiabatic index. We hnd that the g-mode resonances are unimportant for the polytropes 
with N = 0.5 but have some relevant contribution as Fi becomes larger which is mainly 
due to the small values of the overlapping integrals. The inverse is true for the f-modes. The 
f-mode frequency for the case N = 0.5 is smaller and it results that the deviation from the 
point-particle evolution is larger, see Figures 1 to 4. 

In Figures 2 and 3 one can see the differences between the distances of the two stars 
as calculated with radiation reaction -|- tidal interaction on the one side and with radiation 
reaction alone on the other. The differences start oscillating as the orbital motion of the 
two stars approaches frequencies of one half of the oscillation frequencies. The hgures show 
the effects of the /-, ( 71 -, and ( 72 -modes for the binary systems A 0.5 — Gq.s and Ai — Ci. 
Both hgures clearly show that the /-modes dominate the ( 7 -modes. The effects of the tidal 


© 0000 RAS, MNRAS 000 , 000-000 
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effects mildly corrupt the nice sinusoidal form that the standard matched hltering technique 
is supposing. This leads to a smaller correlation and results in an increase of the errors in the 
calculation of the mass parameters of binaries from measured waveforms. Evidently, near 
the high-frequency end of the gravitational-wave-observation window tidal effects should be 
incorporated into the template waveform. This is supported also by the recent works by 
Reisenegger and Goldreich (|1994|) and Shibata (|1994|) . 

For the total phase, see Figure 4, we have obtained a phase difference between the 
(orbital) motion with radiation reaction and the motion with radiation reaction -I- tidal 
interaction of the order ISvr for a binary consisting of two stars of type Ai-Ci (see Table 1) 
and 2 l 7 r for a binary consisting of one type A 0.5 and one type Cq.s star, while in the case of 
a binary with stars of types Aj and Aj the phase difference is quite small and unimportant. 
In both cases these phases have been accumulated during the 100 - 1000 Hz emission stage 
of gravitational radiation while the stars have completed about 300 revolutions around each 
other. 

While the tidal resonances contribute seemingly only during the late stages of the binary- 
systems evolution, the tidal effects accumulate energy on the stellar oscillations during the 
whole period of the coalescence. The transfer of energy becomes signihcant if the system has 
the following characteristics: (i) long coalescing time, which implies small stellar masses since 
the coalescing time is proportional to ~ M~^ for nearly equal-mass binaries, and (ii) small 
ratios rui/ri, because in this case the tidal effects are more signihcant (it results in a larger 
energy transfer from the orbital motion to the oscillations). Even though, the contribution 
of the tidal effects outside the region where resonances get important is in general small. 
For the system Ai-Ci we have mentioned earlier that during the whole evolution of the 
system, from the beginning of the observation window, i.e. lOOHz which correspond to a 
distance of 44. 8 M, the phase difference accumulates to 14.27r, while the evolution of the 
same system from a distance of 20M (frequency of the gravitational waves 335 Hz) to 9.7M 
(end of the observational window, where also the stars touch each other) accumulates a 
phase difference of 12 . 97 r. The corresponding numbers for the system Aq.s-Cq.s are 2 l 7 r and 
18.dvr. This means that 90% of the total phase difference is created during the hnal inspiral. 
This late-stages corruption of the signal due to the tidal effects is also apparent in the work 
of Kochanek (|1992|) where irrotational Roche-Riemann ellipsoids have been used as models 


for the neutron stars. 
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An interesting featnre is the fact that if the two stars are close enough the tidal effects can 
drive them into coalescence even in the case that no radiation reaction is taken into account, 
see Figure 5. The tidally induced plunge region strongly depends on the masses of the two 
stars and on their degree of compactness. The more compact the stars are the smaller their 
separation is for the tidally induced plunge. For the three models that we referred to earlier 
this distance takes the values: 11. 6 M (8.9M) for the system A 2 -C 2 , 16M (8.9M) for the 
system Ai-Ci, 17.4M (8.9M) for the system A 0 . 5 -C 0.5 and 19M (15.9M) for the system Ci- 
Ci (in parentheses are shown the distances for which the stars touch each other). We should 
notice in addition that this tidally induced plunge, in reality, should be more complicated 
since it happens for distances where the two stars have passed already their Roche limit 
( Pildsten fc Cutler 1992|) . Recently Lai et ah (|1993|) , (|1994| ) came to similar conclusions for 
binary systems with also polytropic star components, but they have not used the approach 
which is applied in the present work. The results support the idea of introducing the notion 
of extended bodies in the construction of the waveform templates for the future detection of 
gravitational waves. Nevertheless, we should keep in mind that we have used linear analysis 
and that the “plunge orbit” occurs at small separations at which the linear theory is not 
strictly valid. Thus, although the existence of “plunge orbits” seems to be established the 
numbers that we give might change if non-linear theory is applied. 

Concluding this work we will discuss the effects of our models on the waveform of the 
gravitational waves emitted during the binary coalescence stage. As it can be seen from the 
Figure 4 (remember the end of Sect. 2 ) the waves accumulate very quickly a phase difference 
which seriously corrupts the shape of the Newtonian waveform leading to difficulties in 
the detection procedure if not taken into account in the template waveform (for a numerical 
simulation of the detection procedure and the exact influence on the detection procedure the 
reader is referred to ([Kokkotas et ah 1994a|) ). Although the IPN corrections are weakening 
the high phase differences ([Kokkotas et ah 1994b|) the effects on the phase are very large 
and cannot be overcome. For the evolution of the system Aj-Cj, from 20M to coalescence, 
the phase difference produced by the IPN theory is 327r, while from the tidal effects, as 
mentioned earlier, 2 l 7 r for the system Aq.s-Cq.s and Idvr for the system Ai-Ci. The large 
tidal contributions depend on the choice of the binary system. If we would have chosen 
a system of the type Aj-Aj then the post-Newtonian effects would have been of the same 
order while the tidal ones would have turned out much smaller. This can be considered as a 
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positive result since, in this way, one gets extra information about the compactness of the 
stars. 


The system of equations that have been used in this work would be more precise if 
the radiation loss due to non-radial oscillations of the stars would have been included. 
Nevertheless, to leading order, the above analysis is still valid since the energy stored in the 
stellar pulsations is a fraction of the total energy of the system (notice that the reaction 
Hamiltonian, equation (19), is proportional to the non-spherical-symmetric part of the sum 
of two times the kinetic energy tensor plus the potential energy tensor and, implicitly, to the 
hrst time derivative of that part). Furthermore, the damping time of the stellar oscillations 
which for identical stars is of about the same order as the coalescence time (coalescence 
understood as the stage where the interbody-body separation amounts to some few single¬ 
body diameters only) becomes longer as the stars become less relativistic; this can be proven 
easily from the relations given in Balbinski and Schutz (|1982| ) eq. (1) and/or Reisenegger 


and Goldreich (|1994|) eq. (11). In our work we have taken the most compact star as a point¬ 
like mass while the other star which is less relativistic was the one that has undergone 
oscillations, oscillations with indeed much longer damping time than the coalescence time 
of the binary system. But if the evolution of a binary system is studied from the regime of 
lOHz, since the coalescence time is long, the damping of the stellar oscillations should be 
taken into account. 
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14 Kostas D. Kokkotas and Gerhard Schafer 

TABLES 

Table 1: For i = 2, nine polytropic neutron stars models are listed for three polytropic 
indices N = 0.5, iV = 1 and N = 2. The index on the letters A, B, C shows the polytropic 
index. For these models, which span a wide range of possible neutron stars radii, masses, 
and densities, the /-mode frequency is given in units of Hz. 

Table 2: For i = 2, normalized frequencies for three polytropic-adiabatic index pairs 
are listed together with the corresponding values for the overlap integral Qn 2 - 
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Table 1 


NEUTRON STAR MODELS 

Model 

r (Km) 

m (Km) 

m/M© 

Pc (gr/cm^) 

/-mode (Hz) 

^0.5 

10.105 

2.036 

1.38 

1.17 X 

2212.8 

-Bo.5 

12.533 

1.561 

1.06 

4.69 X 10^^ 

1402.4 

Co.5 

18.162 

1.144 

0.77 

1.13 X 10^^ 

688.4 


10.105 

2.036 

1.38 

2.09 X 10^5 

2607.6 

Bi 

12.533 

1.561 

1.06 

8.40 X 10^^ 

1324.8 

Cl 

18.162 

1.144 

0.77 

2.02 X 10^^ 

809.9 

^2 

10.105 

2.036 

1.38 

7.25 X lO^® 

3742.6 

B2 

12.533 

1.561 

1.06 

2.91 X 10^® 

2372.0 

C2 

18.162 

1.144 

0.77 

7.02 X 10^^ 

1164.3 
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Table 2a 


Eigenfrequencies and Overlap Integrals 

N = 

0.5, ri = 

3.05 

N = 1, 

Ti = 2.05 

Mode 

ujl2[Gm/r^ 

] |Qn2| 

^n2[Gm/r^] 

1 Qn2 1 

P2 

48.2425 

.1155x10“^ 

30.1934 

.2541x10-2 

Pi 

18.8887 

.7475x10-2 

12.3929 

.2588x10-° 

f 

1.0883 

.6238x10+° 

1.5064 

.5580x10+° 

9i 

0.0088 

.3193x10-° 

0.0340 

.1764x10-2 

92 

0.0038 

.5054x10-"^ 

0.0161 

.4250x10-° 

93 

0.0022 

.3106x10-^ 

0.0095 

.1256x10-° 


Table 2b 


Eigenfrequencies and Overlap Integrals 


N 

= 1. ri = 7/3 

N = 2, 

Fi = 5/3 

Mode 

^l2[Gm/r^] 

|Qn2 1 

^l2[Gm/r^] 

|Qn2 1 


61.1329 

.3095x10-2 




40.6315 

.8070x10-2 



P2 

35.1837 

.2204x10-2 

24.0741 

.2287x10-° 

Pi 

14.7537 

.2121x10-° 

11.5549 

.7574x10-° 

f 

1.5112 

.5580x10+° 

3.1133 

.4219x10+° 

9i 

0.1902 

.1110x10-° 

0.5633 

.2105x10-° 

92 

0.0921 

.2607x10-2 

0.2967 

.8596x10-2 

93 

0.0549 

.7634x10-° 

0.1839 

.3844x10-2 


0.1254 

.1821x10-2 




0.0911 

.8946x10-° 




0.0693 

.4494x10-° 
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FIGURES 

Figure 1: The evolution of the separation of a binary system is shown. The system is 
left to evolve from a quite small separation of 2AM or 125Hz orbital frequency. The dash- 
dotted and dashed lines represent the evolution as determined from the Newtonian and IPN 
approximation, respectively, each time including the appropriate radiation reaction terms, 
i.e. the 2.5PN approximation for the Newtonian motion and 3.5PN approximation for the 
IPN motion. The continuous line corresponds to the evolution of the system as dictated 
only from the Newtonian dynamics together with radiation reaction and tidal interaction. 
Slightly inconsistently, 2.5PN (or 3.5PN) means only the radiation reaction part without 
inclusion of the lower ‘integer’ approximations. 

Figure 2: The difference of the separation of the stars is shown as fnnction of time 
considering, additionally to the Keplerian motion, radiation reaction on the one side and 
radiation reaction -|- tidal interaction on the other side. The figure presents the effects of 
the /-, gi- and 5 f 2 -niode of star 1 (model C0.5) with its corresponding frequencies of abont 
688.4, 61.8 and 40.9 Hz. The contribntions of the analogous modes of star 2 (model A0.5) 
are negligible. Time is plotted in terms of R = Rrad- 

Figure 3: The difference of the separation of the stars is shown as fnnction of time 
considering, additionally to the Keplerian motion, radiation reaction on the one side and 
radiation reaction -|- tidal interaction on the other side. The hgure presents the effects of the 
/-, gi-, and 5 f 2 -niode of star 1 (model Ci) with its corresponding freqnencies of 809.9, 121.7 
and 83.8 Hz, respectively. The contribntions from the analogous modes of star 2 (model Ai) 
are negligible. The (m = 2)-resonances for the gi and g 2 occnr near R/M = 39.3 and 50.5, 
respectively. The resonance of the 5 ' 1 -mode can be clearly seen in this Figure. Time is plotted 
in terms of R = Rrad- 

Figure 4: For N = 0.5 and A^ = 1, the phase differences between the evolution predicted 
by Newtonian and IPN approximation are shown. The tidal effects slowly accnmulate a phase 
difference which becomes large only at the final stage of the coalescence leading there to a 
phase difference of a few complete cycles. 

Figure 5: The evolution of the separation of the two stars due to tidal effects only 
(radiation reaction is not included). For each system there is a characteristic distance beyond 
which the system is driven into coalescence without the need of radiation reaction. 
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This figure "fig 1-1.png" is available in "png" format from: 


http://arXiv.org/ps/gr-qc/9502034vl 



Figure 1 
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This figure "fig 1-2.png" is available in "png" format from: 
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This figure "fig 1-3.png" is available in "png" format from: 


http://arXiv.org/ps/gr-qc/9502034vl 
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